Shear-induced vortex decoupling in E^S^CaC^Og crystals 
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Simultaneous transport and magnetization studies in Bi2Sr2CaCu2 0s crystals at elevated cur- 
rents reveal large discrepancies, including finite resistivity at temperatures of 40K below the magnetic 
irreversibility line. This resistivity, measured at the top surface, is non-monotonic in temperature 
and extremely non- linear. The vortex velocity derived from magnetization is six orders of magnitude 
lower than the velocity derived from simultaneous transport measurements. The new findings are 
ascribed to a shear-induced decoupling, in which the pancake vortices flow only in the top few Cu02 
planes, and are decoupled from the pinned vortices in the rest of the crystal. 



Transport measurements are one of the most common 
methods to study vortex dynamics in superconductors. 
The derived resistivity p describes the vortex motion as 
a function of temperature, field, and the applied cur- 
rent. In high-temperature superconductors the situation 
is more complicated due to their high anisotropy and 
layered structure. The corresponding resistivity has two 
main components, the in-plane resistivity p a b and the 
out-of-plane p c , with typical ratio of p c j p a b — 10 4 w the 
normal state of Bi 2 Sr 2 CaCu20 8 (BSCCO) crystals^. As 
a result, the measured resistance R is a non-trivial fuaer 
tion of sample geometry and contact configuration!!!!, 
which is further significantly complicated by the nonlin- 
ear current dependence of p a b and p c . Yet it is generally 
assumed that the physical mechanism that governs the 
dissipation can be described in terms of current densi, 
namely, the local values of p a b and p c are determine* 
by the corresponding in-plane and out-of-plane current 
densities j a b and j c (ignoring possible nonlocal effectsBQ) . 
In this paper we demonstrate that in highly anisotropic 
materials like BSCCO this assumption may not be valid. 
We find that at elevated currents an additional new term, 
the c-axis gradient of the in-plane current dj a b/dz, be- 
comes the dominant parameter in the description of the 
local dissipation. This current gradient induces large ve- 
locity gradients dv a b/dz of the pancake vortices in the 
different Cu02 planes, leading to their decoupling and to 
corresponding dramatic increase in p c and anisotropy 7. 
In the presence of inhomogeneous currents, this mech- 
anism results in fundamental changes in the transport 
behavior, including large measurable reentrant resistance 
well below the magnetically determined irreversibility line 
(IL). 

The studies were carried out on several high quality 
BSCCO crystalsQ, with T c ~ 90 K. Four wires were at- 



tached to the gold pads evaporated on freshly cleaved top 
ab surfaces, as shown schematically in the inset to Fig. 
la. The bottom surface of the crystals, free of electri- 
cal contacts, was attached to an array of 19 2DEG Hall 
sensorsEI, 30 x 30 pm 2 each, allowing simultaneous re- 
sistance R and local magnetization measurements in the 
presence of transport current. In addition, several crys- 
tals were irradiated by very low doses of 5.8 GeV Pb ions 
to produce columnar defects with concentrations corre- 
sponding to matching fields of B<j, = 5, 20, and 60 G. 
In order to focus on bulk vortex dynamics and to avoid 
complications due to surface barriersta, most of the sam- 
ples were prepared in a form of large square platelets 
with electrical contacts positioned far from the edgesEJ. 
Some of the crystals were cut into strip shape for com- 
parison. Although detailed behavior varies from sample 
to sample, the qualitative features described below were 
observed in all the samples regardless of the surface bar- 
riers and the irradiation doses. Here we present data for 
as-grown crystal A of size 1700 x 1300 x 10/ito 3 , which 
was subsequently cut into a strip, and for an irradiated 
crystal B 1000 x 700 x 30pm 3 with = 60 G. 

Our main general observation is displayed in Fig. 1. 
Figure la shows the R(T) of the unirradiated crystal A 
at various applied fields H a \\ c-axis at elevated applied 
current of 30 mA. At lower currents, that are usually 
used in transport studies, R(T) behaves monotonically 
with temperature as shown for comparison by the dashed 
line in Fig. la for 10 mA at 500 Oe. However, when 
I a is increased we find surprisingly a very pronounced 
reentrant behavior, in which R(T) reaches a minimum 
at some characteristic temperature T min , but then in- 
creases again, often by more than an order of magnitude, 
at lower T. At 100 Oe a sharp drop in R{T) is observed 
at the first-order melting transition (FOT), T m , followed 
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FIG. 1. Resistance (log scale) vs. temperature at various 
applied fields, (a) As-grown sample A at 100, 200, 300, 400, 
500, 700, and 1000 Oe; I a = 30 mA. (b) Irradiated sample B 
at 100, 200, 300, 380, 420, 460, 500, 550, 600, 700, 800, 900, 
and 1000 Oe; I a = 15 mA. Inset to (b) shows R(T) in the 
vicinity of T m i„ at 900 Oe on linear scale. 



by a monotonic resistive tail in the quasi-ordered- 
lattice phase, as described previously^. At 200 Oe the 
behavior is still monotonic similar to 100 Oe data. At 
300 Oe, however, R(T) decreases monotonically down 
to our noise level at about 55 K, but then reappears 
again between 46 and 34 K. This reentrant behavior is 
strongly pronounced at 400 Oe and up to fields above 
1000 Oe. In addition, a sharp discontinuity in the deriva- 
tive dR(T) I dT is often observed at T m j n as shown in the 
inset to Fig. lb. The irradiated samples display essen- 
tially similar properties, as shown in Fig. lb, with a few 
consistent differences: (i) The resistive kink at the FOT 
is smeared by the low dose of columnar defects, consis- 
tent with previous studiesEU, (ii) T m j n is generally shifted 
to higher temperatures, and (iii) the anomalous behavior 
appears at lower currents, as shown by the 15 mA data in 
Fig. lb. The general form of R(T) in Fis_l is somewhat 
reminiscent of the "peak effect" observedEj in NbSe2 near 
H C 2, or at the melting transitior£3 in YBa2Cu30y. The 
behavior reported here, however, is much broader and is 
of a fundamentally different nature. 

In the unirradiated crystals the reentrant features usu- 
ally appear above 25 mA, close to our upper bound of 30 
mA, limited by the contact resistance. The enhanced 
pinning in the irradiated samples causes the anomalous 
behavior to begin at I a as low as 10 mA, allowing a more 
systematic study. Figures 2a and 2b thus present R(T) 
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FIG. 2. Resistance at various I a (left axis, log scale) 
and magnetically measured critical current (right axis, lin- 
ear scale, open circles) vs. T for the irradiated sample, H a — 
200 Oe (a) and H a = 500 Oe (b). 



of crystal B at various I a at 200 and 500 Oe respective- 
ly. At low currents R(T) is monotonic and only weakly 
current dependent, and has a rather well-defined temper- 
ature Tr=o & t which R(T) drops below experimental res- 
olution. All these qualities change drastically at higher 
currents. The appearance of non-monotonic R(T) is ac- 
companied by extreme non-linearity. An increase of I a 
by 30% or less may result in enhancement of R by orders 
of magnitude. Furthermore, the same current increase 
causes a decrease of Tr = o by tens of degrees K. 

In order to elucidate the origin of this anomalous be- 
havior we have carried out simultaneous local magnetiza- 
tion measurements in presence of I a as shown in Fig. 3a 
for T = 30 K. The results seem to be paradoxical. Finite 
resistivity should be present only above the magnetically 
measured IL, which in Fig. 3a occurs above 1600 Oe. Be- 
low the IL the resistivity should be immeasurable in stan- 
dard transport measurements. This is indeed the case at 
low I a . However, at elevated currents, substantial R is 
measured concurrently with the hysteretic magnetization 
well below the IL (in Fig. 3a). Figure 3b shows the corre- 
sponding field profile B z (x) obtained by the Hall sensors 
at 400 Oe on increasing and decreasing H a in presence of 
I a . A clear Bean profile is observecO. Fitting-i-his pro- 
file to the theoretical B z {x) in platelet sampled results 
in critical current density J c — 2 x 10 4 A/ cm 2 , which 
translates into total critical current of I c = 4.2 A. Obvi- 
ously, a transport current of 25 mA, which is less than 
1% of I Cl cannot result in measurable resistance in the 
standard critical state models. 
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FIG. 3. (a) Resistance (right axis) and hysteretic magne- 
tization loop in the sample center (left axis) vs. H a at T — 30 
K and I a = 25 mA for crystal B. (b) Profile of magnetic in- 
duction across the sample at 400 Oe on increasing (•) and 
decreasing (O) fields. 

Figure 2 shows I C (T) determined from the Bean pro- 
files together with the R(T) data. The reentrant resis- 
tance always occurs in the region where I a <C I c - We 
can also compare the electrical fields. In the magneti- 
zation measurements of Fig. 3 the slow sweep of H a 
of about 0.3 Oe/sec induces an average electrical field 
of 4 x 10~ 10 V/cm in the sample, which corresponds to 
a very low vortex velocity of 10 -4 cm/sec at 400 Oc, 
for example. The resistively measured electrical field, 
induced by the simultaneous transport current, on the 
other hand, is 1.4 x 10~ 3 V/cm. This translates into 
vortex velocity of 350 cm/sec, which is 6 orders of mag- 
nitude higher. Clearly existence of such a high vortex 
flow rate induced applied by concurrently I a should have 
resulted in completely reversible magnetization, in sharp 
contrast to the data. So paradoxically we have a situa- 
tion in which a transport current, which is two orders of 
magnitude lower than 7 C , results in vortex velocity which 
is 6 orders of magnitude higher than the flux creep vortex 
velocity at I c . 

We emphasize that we do not observe any heating ef- 
fects as confirmed by monitoring the location of the re- 
sistive transitions at T c and T m at various currents. Fur- 
thermore, and more importantly since the two measure- 
ments are carried out simultaneously, the transport and 
the magnetization results should be affected by heating 
to the same extent, and therefore heating effects can be 
ruled out as a possible explanation of the observed large 
discrepancy between the two measurements. 

We describe the observed phenomena in terms of shear- 



induced decoupling and shear-enhanced anisotropy. At 
elevated temperatures in the vortex-liquid phase the in- 
plane current density j a b(z) decreases approximately ex- 
ponentially from the top surface to the bottomu. In sam- 
ples of typical dimensions the current at the bottom is 
about four orders of magnitude lower than at the top, due 
to high anisotropy 7 of BSCCO. As T is decreased the 
pinning becomes more effective. At some temperature 
the current density at the bottom becomes comparable or 
lower than J c . As a result the vortices at the bottom stop 
moving or reduce their velocity substantially, whereas the 
vortices at the top maintain their high velocity since the 
current density there is significantly above J c . As a re- 
sult, the velocity gradient between the planes dv a b/dz is 
increased. The thickness of the vortex pinned layer grows 
as T is decreased resulting in progressively larger veloc- 
ity gradients within the vortex mobile part at the top. 
The enhanced dv a b/dz results, in turn, in shear-induced 
phase slippage between the adjacent Cu02 planes reduc- 
ing the Josephson coupling and leading to decoupling of 
the planestrQ. As a result the p c , and hence the mea- 
sured R, are increased significantly causing the reentrant 
R(T). This process is highly nonlinear, since larger I a 
results in larger dj a b/dz, causing an increase in p c and 7, 
which lead in turn to even shallower current profile and 
thus extremely nonlinear R(I). Since the onset of the en- 
hanced shear is induced by pinning, low dose irradiation 
causes the anomalous behavior to set in at higher T and 
lower I a as seen in Fig. 1. It is important to realize that 
the enhanced p c prevents the lower vortex pinned region 
from effectively shunting the current. In the case of full 
decoupling, for example, p c becomes comparable to the 
normal state p c , which in BSCCO crystals of typical di- 
mensions results in about 10~ 4 f2 resistance between two 
adjacent Cu02 planes. Thus the serial c-axis resistance of 
several decoupled layers in the upper part of the sample 
can be of the order of 10~ 3 £1, thus effectively insulat- 
ing the lower zero-resistance region. Any velocity gradi- 
ent dv a b/dz should, in principle, completely decouple the 
layers, since the time-averaged Josephson phase differ- 
ence reduces to zero. More detailed analysis!!^, however, 
shows that the pancake shear occurs through a 'stick - 
slip' process, which retains a finite Josephson coupling 
at low dvab/dz. Thus the Josephson current between the 
layers can be sustained at low dv a b/dz, and is lost at 
higher rates of pancake shear. 

The described mechanism is conceptually different from 
conventional electrodynamics in which the dissipation is 
determined only by the local vortex velocity, and not 
by the velocity gradient. The underlying physics is also 
very different from the process of nonlocal resistivityeffl, 
in which a force applied to a pancake vortex results in 
a nonlocal drag of pancakes in adjacent layers. This 
drag, which is proportionals to d 2 v a b/dz 2 , induces a more 
uniform vortex flow, in contrast to the described shear- 
induced decoupling which has an opposite effect of am- 
plification of the velocity differences. 

Since in most of the crystal thickness the vortices are 
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FIG. 4. B - T phase diagram of crystal A at 30 mA. T m 
(•) is the first-order melting line. T m in (A) is the position of 
the resistivity minimum, and Tr = o (O) i s the temperature at 
which R falls below the experimental resolution. 

pinned, Bean profiles are obtained in the bulk upon 
sweeping H a (Fig. 3b). At the same time the transport 
current confined to the top planes is sufficient to cre- 
ate vortex flow in these planes, resulting in a measurable 
voltage drop, and explaining the apparent discrepancy 
between the transport and magnetization data. A similar 
mechanism probably occurs when the c-axis properties 
are probed by transverse ac susceptibility and transport 
in the zero-field-cooled state as reported previouslyO. 
Our model also explains the hysteretic R(H) in Fig. 3a. 
The field B in the top planes is determined by the hys- 
teretic Bean profile in the rest of the crystal. As a result, 
the measured R(H a ) displays hysteresis, even though the 
vortices are mobile in the top planes. Consequently, this 
apparent hysteresis practically disappears when R is plot- 
ted vs. the average field B measured by the sensors, 
rather than vs. H a . 

Figure 4 summarizes the B — T diagram of sample A 
at 30 mA. A similar diagram is obtained for the irradi- 
ated crystal B. The resistance vanishes within our reso- 
lution at Tfl = 0j the position of which is highly current- 
dependent. The resistively determined FOT extends to 
fields of about 500 Oe in this sample. The anomalous 
reentrant R(T) is observed in the broad region between 
Tr = q and T m j n , where vortices are pinned in the bulk. It 
is interesting to note that the decoupling process occurs 
first in the entangled vortex-solid phase above 500 Oe, 
where vortex cutting is apparently relatively easy. As I a 
is increased, the shear-induced decoupling expands also 
into the quasi-ordered-lattice phase below T m as is the 
case in Fig. 4. The non-linearities expand also partially 
into the liquid phase above T m i n and T m , where the shear 
process can still modify p c at higher currents. 

One can evaluate the number of the Cu02 planes par- 
ticipating in the flow of vortex pancakes. In Fig. 2b at 30 
K, for example, the total critical current determined from 
the Bean profiles is I c = 3.7 A. This translates into crit- 
ical current of about 0.18 mA per (double) Cu02 plane of 
the sample. Thus, the top 140 planes would be sufficient 
to carry the entire / a =25mA without observable dissipa- 
tion. In order to obtain the large measured resistance the 



actual current in the planes has to be significantly above 
0.18 mA. Furthermore, the lower part of the sample also 
carries a substantial portion of I a . We therefore con- 
clude that only a few tens of Cu02 planes at the top of 
the crystal participate in the described process. Current- 
induced decoupling was recently .analyzed theoretically 
in absence of Josephson couplings! It was found that 
above some critical value of the current applied to the 
top plane, slippage between magnetically coupled planes 
occurs, leading to an abrupt increase in the vortex ve- 
locity in the top plane. Since in this model p c is in- 
finite, there is no corresponding current redistribution. 
Computer simulations!^ of Josephson-coupled planes also 
show sharp current-induced decoupling. Vortex cutting 
was also observed in YBCO crystals^ elevated currents 
in flux transformer measurements! 1 9 rr°l. The much lower 
anisotropy of YBCO, however, causes only a small change 
in the transport properties. In BSCCO in contrast the 
shear-induced decoupling leads to variations by orders of 
magnitude in the observed behavior. 

In conclusion, we reveal large apparent discrepan- 
cies between simultaneously measured magnetization and 
transport properties, which reveal an important feature 
of vortex dynamics in the presence of inhomogeneous cur- 
rents. In layered superconductors at elevated currents 
the c-axis resistivity becomes a strong function of the lo- 
cal current gradient dj a b/dz, instead of being determined 
just by the current density. This mechanism results in a 
highly nonlinear reentrant resistance well below the irre- 
versibility line. 
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